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Abstract 

We derive a continuum field theory for the Majumdar-Ghosh model in the large-S" 
limit, where the field takes values in the manifold of the SO{3) group. No topological 
term is induced in the action and the cases for integer spin and half-integer spin appear 
to be indistinguishable. A one-loop /3— function calculation indicates that the theory flows 
towards a strong coupling (disordered) phase at long distances. This is verified in the 
large- A^ limit, where all excitations are shown to be massive. 
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1. Introduction 

The study of quantum spin models - quantum antiferromagnets in particular - has 
attracted a great deal of interest [ 1 - 4 ] , ever since it was recognised that the ground 
state properties of the Heisenberg antiferromagnet (AFM) and its generalisations in two 
dimensions could be relevant to high T^. superconductors [ 5 ]. One widely used method 
of analysis which is applicable in the long-wavelength limit has been the continuum field 
theory approach, wherein an approximate field theory is derived for the quantum AFM in 
the \aige-S limit and the properties of the ground state are obtained using standard field 
theory techniques. 

For one- dimensional quantum spin chains, the field theory approach led to some in- 
teresting and unexpected results. By mapping the Heisenberg AFM to a non-linear sigma 
model defined on the manifold of S"^ , Haldane [1,2] argued that an integer spin system 
would have a ground state with massive excitations and exponential disorder, but a half- 
integer spin system would have massless excitations and algebraic order. The distinction 
between integer spins and half-integer spins was shown to be caused by the presence of a 
topological term in the action. In two dimensions, however, for the Heisenberg AFM on 
bipartite lattices, the current understanding [ 3 , 4 ] is that no topological term is induced 
in the effective long-distance action and the system always exhibits long-range Neel order 
(at zero temperature) with massless excitations. 

However, for frustrated spin systems in both one and two dimensions, the ground states 
have not yet been conclusively determined. It has been plausibly argued that the frustrated 
J — J' model (with both nearest-neighbour and next-nearest-neighbour couplings) on 
bipartite two-dimensional lattices would have a fiux- phase ground state [ 6 ]. Another 
example of a two-dimensional frustrated model - the Heisenberg AFM on a triangular 
lattice - has been studied [ 7 ] in the field theory limit. It was shown there that no 
topological term is induced and no distinction between integer and half-integer spins is 
observed. 

In one dimension, the Majumdar-Ghosh model [ 8 ], which has both nearest-neighbour 



and next-nearest-neighbour couplings with a particular ratio of their strengths, is a specific 
example of a frustrated spin system. This model has been solved for S = 1/2 and has 
been shown to have doubly degenerate valence bond ground states with massive excitations 
[ 8 ]. However, frustrated spin models in one dimension have not yet been studied from 
a field theoretic point of view. In this paper, we shall address ourselves to a study of 
the Majumdar-Ghosh model in the large-S", long-wavelength limit. We shall derive an 
appropriate field theory for the model and analyse it. In Sec. 2, we briefiy review the results 
from the spin wave, field theory and large-A^ approaches to the unfrustrated Heisenberg 
AFM. In Sec. 3, we study the Majumdar-Ghosh model in detail. We derive a spin wave 
theory about the classical ground state using Villain's approach [ 9 ] and show that the 
model has three spin wave modes, with two of them having the same velocity and the 
third having a higher one. We derive a continuum field theory for the model in the large- 
S, long-wavelength limit and discuss its symmetries. No topological term is induced in 
the action. Hence the long-wavelength theory appears to be indifferent to the distinction 
between integer and half-integer spins. We then perform a one-loop /9— function analysis 
of the field theory which indicates that the system fiows towards a disordered phase at 
long distances. Finally, we generalise the field theory to an appropriate large-A^ field 
theory and show that a saddle-point approximation indicates that all the excitations are 
massive. In Sec. 4, we summarise the evidence offered in the earlier sections and present 
our conclusions. 



2. The Heisenberg Antiferromagnet 

In this section, we shaU briefly review [1,2] the study of the simplest spin chain - 
the nearest-neighbour Heisenberg AFM - given by the Hamihonian 

H= J Y^ S,-Sj (2.1) 

<i,j>&n.n 

with J > 0. The classical ground state (Neel state) of this model has neighbouring spins 
antiparallel. The standard spin-wave analysis can be carried out about this ground state 
using the Holstein-Primakoff transformation [ 10 ] 

SI = S - a\ai 

S- = a\ {2S - a\a,f''' ^ V2S aj (2.2) 

S+ = {2S - aJa,)^/2 ^. ^ ^/^ ^^ 

where the second approximate equality is true in the large-S" limit. Two spin-wave modes 
are found at k = and k = n satisfying the relativistic dispersion E = ck with c = 2JSa. 
However, it is well-known [ 11 ] that there exists no long-range order in one dimension, 
due to infra-red divergences in the theory. Hence, a small fluctuation analysis about an 
ordered state is plagued with inflnities and breaks down in the disordered phase. 

The next approach is the continuum fleld theory approach initiated by Haldane [ 1 ], 
which led to the startling result that integer spins have massive modes and exponential 
disorder, whereas half-integer spins have massless modes and algebraic order. Here, how- 
ever, we shall follow Affleck's [ 2 ] method to obtain the fleld theory. Continuum flelds (p 
and 1 are deflned as 



and (p2i+i/2 = ^ = nx) 

where a is the lattice spacing. These variables satisfy the identities 

1-^=0 



(2.3) 



and (j) =1 + - — 



so that in the large- S" Umit, (jP = 1, - i.e., (p takes values on the manifold of the sphere 
5"^. The commutation relations for the spins imply that the fields satisfy the commutation 

relations given by 

[ /„, //3 ] = «ea/37 ^7 S{x -y) 

[ /a, (t)p ] = ie^pj (pj S{x - y) (2.5) 

[ (/>a, (?^/3 ] = 

which indicates that \{x) can be identified as the angular momentum of the field (f) con- 
structed as 1 = (/) X yf, (where yf is the canonically conjugate momentum to (j)). The 
Hamiltonian in Eq. (2.1) can be rewritten in terms of these fields and expanded in a 
derivative expansion. In the long- wavelength limit, we only need to keep terms upto two 
space-time derivatives. In fact, for 1, since it is already proportional to tt, we only need to 
keep terms upto 1^, 1 and 1', where the dot and prime denote time and space derivatives 
respectively. In this limit, the Hamiltonian can be reexpressed as 

H = jdxi'^il+^^'f + ^<?''] (2.6) 

where c = 2JSa is the spin-wave velocity and the coupling constant g^ = 2/S, so that 
large-S" corresponds to weak coupling. 

This Hamiltonian in Eq. (2.6) follows from the Lagrangian given by 

where the field (p is subject to the constraint (p'^ = 1. Setting the spin- wave velocity c = 1, 
the action can be written in a relativistically invariant way as 

S = jd^x[^dj-d^$ + l-e^-$.djxdj] (2.8) 

with 9 = 2nS. The 9 term can actually be shown to be a total derivative, which affects 
neither the classical equations of motion nor the perturbative Feynman rules. However, it 
does have a topological significance. On compactified Euclidean space, both x^ and (p can 
be represented as points on a sphere S'^. The functional 

Qi^) = ^ fd^^ e^"8 [f djx dj] = Q (2.9) 



measures the winding of the sphere S'^{(j)) onto the sphere S'^{x^) and is always quantised 
as an integer for any field configuration, since 112(5'^) = Z. Thus, for integer values of 
5", the contribution of the topological term to the path integral is given by e"^'^* '^ = 1 - 
i.e., it is irrelevant. The action in Eq. (2.8) without the topological term is well-known to 
describe a theory which flows to a strong coupling, massive phase with masses of the order 
e~'^^ . But for half-integer spins, the topological term e"^'^*'^'^ = e~™^ is either +1 or —1 
depending on the value of Q and is certainly relevant. In that case, the topological term 
cannot be ignored. The current understanding [ 2 ] is that the action in Eq. (2.8) with 
a topological term is actually massless and describes a conformally invariant field theory 
with algebraic order - i.e., with a power-law fall-off of correlation functions. 

The result that the field theory in Eq. (2.8) describes a massive phase can also be 
confirmed by generalising the S"^ manifold of the (j) field to an S^ manifold and then 
studying the large- A^ limit of the field theory. The appropriate Lagrangian is given by 

C= ^[dj-d^$ + tX{$^ - 1) ] (2.10) 

where the constraint has been enforced using an explicit Lagrange multiplier field. We can 
now integrate out (p and obtain 






iX 
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5eff(A) = — - I (fx — + trln(-a2 + iA) . (2.11) 



1 


1 , A 




= In 


9^ 


27r m 



In the large- A^ limit, we can ignore fiuctuations of A and evaluate iSeff (A) at the saddle-point 
defined by dSes/dX = 0. We find that the saddle-point equation gives 

(2.12) 

where A is an ultra-violet cut-off. This determines the mass parameter m in terms of the 
cut-off and the bare parameter g as 

m = Ae-^'^/^' = Ae"'^^ (2.13) 

which agrees with the expectation from the one-loop renormalisation group flow. 



Thus, the current understanding of the Heisenberg AFM spin chain is that both for 
integer and half-integer spins, there exists short-range order, but for integer spins, at large 
distances, long-wavelength fluctuations destroy the order. For half-integer spins, even in 
the long-distance regime, excitations are massless and there exists algebraic order. 



3. The Majumdar-Ghosh Model 

The Hamihonian for the Majumdar-Ghosh model is given by 



H^ J J2^^- [S^+i + ^S,+2] (3.1) 



2 

The classical ( 5* ^ oo ) ground state must have any three neighbouring spins adding upto 
zero. This can be seen by recasting the Hamiltonian in Eq. (3.1) as 

if = - ^ ( S, + Si+i + Si+2 f + constant (3.2) 

i 

For this to hold all along the chain, the spins must lie on a plane and must successively 
twist by 27r/3 as shown in Fig. 1. 

For S = 1/2 , this model has dimerised ground states which are just products of local 
singlet pairs. Since each spin can form a singlet by pairing with its neighbour to the left 
or with its neighbour to the right, it is clear that two such ground states are possible. 
These states were proven to be exact ground states by Majumdar and Ghosh. Later, it 
was shown [ 12 ] that these two states are the only ground states and that there exists a 
gap in the spectrum. 

However, the model remains unsolved for all higher spins. Here we shall obtain the 
ground state properties of the model in the large-S", long-wavelength limit. Even though 
naive spin- wave analysis is inadequate because of infra-red divergences in one-dimensional 
systems, we perform a spin wave analysis to obtain the spin wave modes and their velocities. 
The purpose of this calculation is to have a check on the field theory formulation. This is 
done in Sec. (A). In Sec. (B), we derive the field theory, discuss its symmetries and show 
that it reproduces the spin- wave calculations in the appropriate limit. In Sec. (C), we 
study the renormalisation group flows of the coupling constants. The global symmetries 
of the fleld theory allow for four independent coupling constants. We study the flow of 
these couplings by deriving the one-loop /?— functions and integrating them numerically. 
We obtain the scale at which a transition to the strong coupling (disordered) phase occurs. 
Finally, in Sec. (D), we show that in the large- A^ limit, all excitations become massive. 



This confirms the fact that the theory is in the disordered phase in the long-wavelength 

limit. 

Sec. (A): Spin- wave Analysis 

Since the classical ground state of the Majumdar-Ghosh Hamiltonian is planar, it is 
more convenient to obtain the spin-wave theory by parametrising [ 13 ] the spin operators 
in terms of Villain's variables [ 9 ], rather than the usual Holstein-Primakoff variables. The 
spin variables at each site are expressed in terms of two conjugate variables Sf and (pi as 



sr = 



(s + lf - (st + ^f 



1/2 



{s + l^ - {st + ^f 



1/2 



e 



-t(Pi 



(3.3) 



Sf-' ^ 



i d(fi 
The periodic operator ipi { ipi = ipi + 27r ) satisfies the commutation relation 

[ ^. , ^ ] = ;^ '^^^ ■ (3.4) 

In the large-S" limit, the discreteness of Sf / S and the periodicity of ipi can be neglected 
and SI and ipi can be thought of as canonically conjugate, continuous momentum and 
position operators. In the classical ( 5" ^ oo ) limit, in fact, Sf / S ^> Q and (fi is fixed 
at its classical value ipi , which is the angle made by the classical spin relative to the x- 
axis. The large- 5" approximation involves a systematic 1/ S expansion around this classical 
ground state. 

For the Majumdar-Ghosh Hamiltonian, from Fig. 1, it is clear that (pi — (pj for any 
two neighbouring spins is always 27r/3 . Let us now compute the spin wave spectrum for 
the Hamiltonian by expanding the square roots in Eq. (3.3) and keeping terms only upto 
order 1/ S , and by expanding the angles (fi to quadratic order about pi - i.e. we have 



St - e^"^ {S +]-) {l + ^9. " | ) (l " ^^^ ^' 



2 ' ' 2 ' \ 2 (5+i )2 

iw. .. ^^A iiS! + l)^ 



(3.5) 



s- ~ e-^'^^ {S + -){l-ie, - ^ ) 1 



2 ' ' ' 2 ' V 2 (5+i )2 



where 6i describes in-plane fluctuations (fluctuations about the angle ipi ) and S^ / S 
describes out-of-plane fluctuations. Substituting these operators in the Hamiltonian in 
Eq. (3.1), we flnd that 

H = jY,[st su, - sf cos e - ^ {e, - e, f cos e ] 

2 (3-6) 

i 

where Q = (pi — (pj = 27r/3 . Upon Fourier transforming, we can write the Hamiltonian in 
the standard oscillator form as 

H = J[Y. %^ + IkO.O.,] (3.7) 



2 m„ 2 

q ^ 



where 



u{q) = J-^ = - J^cos e J (3-7J (3-7g/cose ) (3.^ 



rriq 



with 

jq — 2 COS q + cos 2g. (3.9) 

Here, — tt < g < tt . Also note that uj{q) = uj{—q) . The spin-wave spectrum is now 
obtained by expanding uj{q) about its zeroes -i.e., about any qo where uj{qo) = - and 
identifying its velocity as du/dq \q=q^^ ■ It is clear from Eq. (3.8) that uj{q) has three 
zeroes at g = 0, 27r/3 and — 27r/3 within the flrst Brillouin zone. At g = , the 
spin- wave velocity is given by 

c{q = Q) = J S a W ^ (3.10) 

and at both q = 27r/3 and q = — 27r/3 , it is given by 

/27 
c (g = ± 27r/3 ) = J S a \ —. (3.11) 

V 8 

Hence, we have three spin-wave modes, two with the same velocity and one with a higher 
velocity. 

10 



As was mentioned in Sec. 2, spin- wave theory about an ordered state is not appro- 
priate in one dimension, because infra-red divergences in the theory drive any ordered 
state towards disorder. However, the spin-wave velocities computed here will serve as 
a consistency check on the field theory of the model, which will be derived in the next 
subsection. 

Sec. (B): The Continuum Field Theory 

The first task in the derivation of a continuum field theory for the Majumdar-Ghosh 
model is the identification of 'small' and 'large' variables from a local set of spins, analogous 
to the 1 and (j) fields in Eq. (2.3). Since, the classical ground state is a three sub-lattice 
Neel state, we define field variables involving three local spins as 

1 _ ( Ssi-i + Ssi + Ssi-i-i ) 

hi — 7. 

6 a 
(M,. = '^"'"^/"^'' (3.12) 

.7* N ( Ssi-i + Ssi+i — 2 Ssi ) 

[(p2)3i — ^-^ 

where a is the lattice spacing. Note that we have grouped the spins as (3z — 1, 3i, 3i + 1). 
The two other possible groupings are (3z — 2, 3i — 1, 3i) and (3i, 3i -|- 1, 3i + 2). All three 
groupings lead to the same field theory. (We shall verify this later and see how these three 
possibilities are related to symmetries in the field theory.) 

For the classical configuration in Fig. 1, we see that 1 = and the fields (pi and (p2 
satisfy the conditions given by 

<^ = ^ = 1 and (^1 . (^2 = 0. (3.13) 

We may therefore expect 1 to be a 'small' variable (proportional to a space-time derivative 
in the long-distance limit) and (pi and (p2 to be the two 'large' and non-linear variables, 
describing three degrees of freedom (since (/>i ■ (/>2 = 0). The actual identities satisfied by 

11 



these fields are given by 



^? + </>! 



2 + 



2a^l 



2i2 



^2 



>1 ■<P2 



s S^ 

^</>l • 1, 



(3.14) 



which reduces to Eq. (3.13) in the long-distance (where | a\ \ « | (/>i |, | (/>2 | ) and 
large-S" hmit. Moreover, the 1, (f)i and (f)2 fields satisfy the commutation relations given by 

[ la (x), li3 {y)] = i ea/37 l-yi^x) d{x - y) 

[ la (x), (pap (y)] = i ea/37 (Pa-yi^) S{x - y) (3.15) 

[ (t)aa ix),(pbf3 (y)] =0 

where a, b take the values 1 and 2. These relations can be easily derived from the com- 
mutation relations of the spin operators by replacing the lattice 5-functions by continuum 
ones - i.e., dij/3a — > d{x — y). Thus, 1 can be identified as the angular momentum 
operators of the fields (/)i and (/)2- 

We can now derive the continuum Hamiltonian by rewriting Eq. (3.1) in terms of the 
fields 1 , (j)i and (/)2 and then Taylor expanding the fields upto second derivatives of (/)i and 
4>2 and first derivatives of 1, and by replacing the sum over sites by an integral over space. 
Thus, in terms of the fields, we obtain 



H 



dx 



2 ^ ^ ' 



where 



+ v^^^' 



27 



" + $r 



(3.16) 



and g^ = v6/S . 
As before, we see that large-S" corresponds to the weak coupling, perturbative regime. 
Since analysis of the field theory is simpler from a Lagrangian formulation, let us deduce 
a Lagrangian for this field theory. We define a third unit vector 



>3 = <Pl X <P2 

12 




with 

<?i-^3 = '?2-<^3=0 , ^'1 = 1. (3.17) 

This field satisfies the commutation relations expressed in Eq. (3.15) for a, 6 = 3 as well. 
Let us now define an orthogonal matrix R whose entries are 

(</'ll 4>21 4>31 \ 
'/>12 fe <^23 • (3.18) 

(pSl </>32 (pSS J 

Thus, the three columns of ^ are given by the vectors (pi, (j)2 and (j)^ respectively. Notice 
that since R is orthogonal { R R = RR =/), it only has three independent degrees 
of freedom. If we now introduce a diagonal matrix 

(3.19) 

it is clear that the second term in the Hamiltonian in Eq. (3.16) is proportional to 
tr (^ Rl2)- But the problem is to figure out a time derivative term in the Lagrangian 
involving the R fields that would lead to the first term in the same Hamiltonian. 

To figure this out, we note that our problem involves an 5'0(3)-valued field at each 
space-time point. Hence, we need the method of quantisation on a group manifold which 
is derived in Ref. [ 14 ]. Following their work, we can show that the required Lagrangian 
is given by 

£ = -^ tr {ff R) - -%tr {g^ g h) + ^ ^i ■ ^'i x c^^. (3.20) 

Acg-^ 2 g^ v^3 

Let us verify that the Hamiltonian in Eq. (3.16) can be derived from this Lagrangian [ 14 ]. 

R may be parametrised by three numbers ^a, which are local coordinates on the manifold 

of 5*0(3). Thus, the derivatives of ^ are 

d^,R = a^^a d^R (3.21) 

where da = d/da ■ The momenta tTq canonical to the ^a can be found from the Lagrangian 
as 



_ d C 1 _. „T. ^ ^T _ „_ S 

'a — 



7r„ = —^ = —^ tr [ (a, R^) R + g (dc^R)] + ^ (</>! X (j)[) ■ da(j)i (3.22) 
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leading to the Hamiltonian 

H = f dx [-^ tr {R^ R) - -^tr {g^ g h) ]. (3.23) 

Now, the tTq are locaUy defined momenta conjugate to the locally defined coordinates ^a • 
The globally defined momenta (conjugate to the globally defined coordinate variables R ) 
are given by 

L = T^p Np^ (3.24) 

where {dp R ) Npa. = i T^ R and T^ are the generators of SO (3) , so that 
{Ta)p^ = i Coip-y ■ Hence, we find that 

^« = T7^*^ (^^ T^R) + ^ (<?! X (^/) ■ (T, $i). (3.25) 

^cg V3 

Here, we have used the fact that the first column of the matrix Tq, R is given by Tq, (pi . 
The construction of /q, guarantees that 

[ /a (x), Ip (y)] = i ec,p^ lj{x) 5{x - y) 

(3.26) 
[ /a (x), R{y)] = Tc, R (x) 5{x - y) 

which is identical (in component notation) to Eq. (3.15). Moreover, by substituting for R 
in Eq. (3.23) in terms of l^ using Eq. (3.25), we can verify that the Hamiltonian in Eq. 
(3.23) is identical to that in Eq. (3.16). 

Let us now study the Lagrangian given in Eq. (3.20), which we have just proved is the 
continuum field theory of the Majumdar-Ghosh model in the large-S" limit. Notice that 
the last term in Eq. (3.20) is a total derivative, and therefore, has no effect perturbatively. 
In fact, we shall now argue that this topological term actually vanishes for all smooth 
field configurations in Euclidean space, and has no non-perturbative effect either. To 
demonstrate this, we observe that the Euclidean action (setting c = 1) is given by 

d'x {—^ (</>!+ </>2 + </'3 ) + ^ (</>/ ' + </'2' ' )} - ^ Q [</>!], (3.27) 

where Q \(p]\ is defined as Q [^] was defined in Eq. (2.9). Finiteness of 5"^ implies that (/>!, 
(p2 and (/)3 must go to constants at infinity, so that space-time is compactified to 5"^ . The 

14 



functional Q [cpi] must still be an integer since (pi also lies on 5"^ {(p'^ = 1 ). However, 
since (/)i is actually part of an SO{3) matrix R and n2(5'0(3)) = , Q[4>i] must actually 
vanish. Let us demonstrate this explicitly. We define a vector 

(P (0) = (p^ COS e + (p2 sin e. (3.28) 

Clearly, Q [(j){0)] is also an integer, since (j){9) also takes values on 5"^ . However, since 
(^(tt) = — ^(0) , Q[(j){w)] = — Q[(/){0)] . Since Q[(/)] is restricted to integer values, it cannot 
change when 9 is continuously varied. Hence, the only consistent possibility is that Q[4>] 
must vanish for all 9. 

Thus, we have the important result that there is no topological term in the field theory 
of the Majumdar-Ghosh model in 1 + 1 dimensions. Notice that the proof of the vanishing 
of the topological term needed the global existence of two orthonormal vectors (pi and 
(f)2 5 whereas in Sec. 2, the field theory of the Heisenberg AFM involved only one globally 
defined field (p . Hence, the field theory defined in Sec. 2 did possess a topological term 
which distinguished between integer and half-integer spins. Here, however, we can work 
with the simpler 5*0 (3) -valued field theoretic action defined by 

Se= f d^x[^ tr (R^R) + ^ tr {R '^ g h) ] . (3.29) 

This action is invariant under the global symmetry R — * ARB , where A is an 5*0(3) 
matrix and S is a block diagonal matrix of the form 

(cos 9 sin 9 0\ 
-sin 9 cos 9 (3.30) 

1/ 

Thus, the theory has an SO{3)l x SO{2)ii symmetry, where the 'L' index mixes the rows 
and the 'R' index mixes the columns. The 5(9 (3) l denotes the original spin symmetry of 
the Hamiltonian generated by the angular momentum \{x). Equivalently, we may think of 
this symmetry as the symmetry of the 'top' R about a 'space-fixed' set of axes that affects 
the three vectors (pi , (j)2 and (p^ in the same way without mixing them. On the other 
hand, the SO{2)fi is a rotation about the 'body- fixed' axis ^3 which mixes (pi and (p2 ■ 

15 



This symmetry is not present in the original spin Hamiltonian in Eq. (3.1) but only arises 
in the long- wavelength field theory developed around the ground state of Fig.l. We now 
recall the statement following Eq. (3.12). To define the fields, the spins can be grouped in 
threes in three different ways. These three ways are related to each other by rotations in 
the ((/>!, (/)2 ) plane by the angles 0, 27r/3 and 47r/3 . So we should certainly have expected 
a cyclic Z3 symmetry in the field theory. However, the vanishing of the topological term 
in the action appears to have allowed an enlargement of the symmetry to a full continuous 
SO [2) symmetry. 

Finally, let us derive the spin-wave spectrum from the Lagrangian and check whether 
we get the right spin- wave velocities. The classical ground state is clearly given by the con- 
figuration S3. = 5 (0, 1, 0), S3.-1 = 5 (-^3/2, -1/2, 0) and 83^+1 = ^ (v^/2, -1/2, 0), 
which can be expressed in terms of the matrix R by the identity matrix i.e. R = I. Small 
fiuctuations about this ordered state can be parametrised by 



^=7 1 -a \ (3.31) 




where a , (3 and 7 are small fluctuation fields. In terms of these fields, the Lagrangian is 
given by 



': = ^ (4' + -3' + i') - ri ("" + f>'^ + 27" ) . (3.32) 



Thus, the a and (3 modes have the velocity c and the 7-mode has the velocity \/2c in 
agreement with the spin-wave velocities in Eq. (3.10) and (3.11) since c = JSa (27/8)^'^ 
from Eq. (3.16). From the matrix R in Eq. (3.31), it is clear that the 7-mode describes 
fiuctuations within the plane of the classical vectors (j)i and (j)2 - i.e., within the plane of 
the classical spins. The a and (3 modes describe out-of-plane fiuctuations. 

Thus, we have proved that the field theory of the Majumdar-Ghosh Hamiltonian 
reproduces the standard spin- wave results in the weak coupling limit. 
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Sec. (C): One-loop /3- function and the flow of coupling constants 

In the previous section, we derived an effective low-energy field theory for the 
Majumdar- Ghosh model and showed that it reproduced the standard spin-wave scenario 
by expanding about the classical, ordered ground state. However, more generally, we can 
study the long-distance properties of such a field theory by applying the renormalisation 
group. The SO{3)l x SO{2)fi global symmetry of the action in Eq. (3.29) implies that at 
any scale, the effective (Euclidean) Lagrangian can have four different coupling constants 
and can be written as 

^i^ = ( ^ - A )tr (^^ ^) + ( ^ - ^ )tr {R^ R h) 

At microscopic distances of the order of the lattice spacing a, where we derived the field 
theory, we have 

2 2 2 o 2 2 

9l =92 =93 =2^4 =9 

with 

g^ = y^/S . (3.34) 

(Compare Eq. (3.29) and Eq. (3.33) ). But these values change as we move to larger 
distance scales in accordance with the renormalisation group equations. (Note that an 
SO{3)l X SO{3)ii symmetric Lagrangian would have had gi = g2 and gs = g4 )■ 
The unusual parametrisation in Eq. (3.33) has been chosen for convenience in studying 
the evolution of the small fluctuations a, P and 7, so that the spin-wave (Minkowski) 
Lagrangian takes the simple form given by 



£ = -L (d2 + /32 ) - -i^ (a'2 + /3'2 ) ^ 1 ^2 _ 1 y. _ (335) 



^g-f -' ^ ^ ' 2gi ^" + f' ) + 2gl ^ 2gl 



From this Lagrangian, we read off the velocity of the a and P modes to be gi /gs and 
that of the 7-mode to be (72 /g4 . 

To study the flow of the coupling constants to the long distance regime, we calculate 
the four /3-functions 

P{9f ) = ^ 9ny) (3.36) 
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where y = ln(L/a ) is a measure of the distance scale L. This is done using the background 
field formalism [ 15 ]. We expand the field R{x) as 

R{x) = Roix) e*''"(")'^« (3.37) 

where ^o("^) ^^ ^ slowly varying field, which we take to be a solution of the Euler- 
Lagrange equations of motion of the Lagrangian in Eq. (3.33), and rjoiix) are rapidly 
varying fields. We then integrate over the r]a fields to obtain an effective action for 
^0 7 from which the /9-functions are obtained. (The field R{x) has to be expanded as 
R{x) = Rq{x) exp{irjoi{x) T^) rather than R{x) = exp{irjct{x) T^) Ro{x) in order to 
maintain the SO{3)l x SO{2)fi symmetry for the effective action for ^q ). 

On expanding the right hand side of Eq. (3.33) to second order in rjcdx) (the first order 
terms vanish because ^q extremises the action), we find that the Euclidean propagator 
MQ,^(k) = (?7Q,(k)?7^(— k)) in momentum space k = (/cq , ki ) has a diagonal form with 



gi ki + gf ki 
and M33 = 2 l2 !i "^ 2 l2 

94: '^0 + 92 '^l 



(3.38) 



In addition, there are two kinds of vertices between the ^g'^^lds (denoted by dotted lines) 
and the r^-fields (denoted by solid lines) as shown in Figs. 2 (a) and 2 (b). The vertex 
in Fig. 2 (a) has one derivative of ^q and one derivative of 77. (Since the momentum 
flowing along the Rq line is negligible compared to the momentum flowing along the tj 
lines, we have shown the ry-momenta to be k and — k respectively). The Feynman rule for 
this vertex is given by 



ra/3 = -ike 



—iki 



1 1 • T 

9i '^92 

1 1 • T 

+ (— — — ) tr ^Q ^Q ( Tq /2 T/3 — Tp I2 Toi) 
92 9i 

(~ ~ TTT ) *^ ^0 Ro {TaTp - Tf3 Ta ) 

93 ^y4 
11 ,j^ 

+ (— — — ) tr ^Q Rq { Ta I2 Tp — Tp I2 Tq, 
9^ 93 
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(3.39) 



The vertex in Fig. 2(b) has two derivatives of ^q ^^^ '^o derivatives of r], and its Feynman 
rule is given by 

1 1 • T ■ 1 1 

Ta/? = -(—-—) tr {^0 ^0 (Ta I2T/3 - - Ta Tf3 h - -h T^, Tp )} 

- (^ - 4 ) tr {gl g^ (T« I^Tp - \ T, Tp h - \l2 T„ Tp )} . 
94 93 II 

The one-loop contribution to the effective action is therefore obtained by evaluating 
the diagrams in Figs. 3 (a) and 3 (b), which arise from the vertices in Figs. 2 (a) and 
2 (b) respectively. Both the diagrams in Fig. 3 have symmetry factors of 1/2 . The 
momentum integrals are cut off at short distances by the lattice spacing a and at long 
distances by the length scale L at which we wish to evaluate the effective action - i.e., 
-^"^ ^ ^0? ^1 ^ a~^ . It is also convenient to use identities like 

( <?! ■ 4 )' = ^ ( <?? + <?2 - ^3 ) (3.41) 

which follow from the constraints in Eq. (3.13) and Eq. (3.17). We then obtain the 
/9-functions given by 



ii{9\ ) = 


9\ 

27r 


r 9i 93 94 2 1 1 

[ 2 / , -. \ 9i 93 [ 2 2 
9i [9i 94 + 92 93) 91 9i 


ii{9l ) = 


9\ 

27r 


r 3 / 2 1,2 /I 1 N 1 

[9i93[ 2 2 ) + '^9i93[ 2 2 ) ] 
9i 92 92 9i 


/5(^i ) = 


9\ 

27r 


r 9I 91 92 2 1 1 

2 / , \ + 9i 93 [ 2 2 
91 [91 94 +92 93) 93 91 


(i{9l ) = 


9\ 

27r 


r 3 . 2 1 ,2 .1 1 ^ n 

[^3^1 (9 9 ) + 2(71 (73 (2 ^2 )]• 

93 94: 94 93 



(3.42) 



Let us now integrate the /3-functions numerically starting from their values at y = 
given in Eq. (3.34). The one-loop /3- functions have the property that they are invariant 
under the rescaling 

y — > \y and gf — > gf/\ (3.43) 

which implies that the /3-functions for different values of g'f (or equivalently 1/5" ) are 
related by a scaling law. This also means that if one or more of the coupling constants go 
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to infinity at any length scale y = yo , then the product yo g^ (where g^ = V^/S ), must 
be a number independent of the value of g^. 

From our numerical study of the renormalisation group fiows, we find that the coupling 
constants run in such a way that the velocity of the 7-mode increases and the velocity of 
the a and f3 modes decreases. (This behaviour was verified over a large range of initial 
values starting from g'^ = 0.001 through g^ = 1.0 ). Although all four couplings increase 
monotonically, the coupling (72 goes to infinity first, and this happens at a length scale 
given by 

yo / = 4.5 (3.44) 

(At that point, the other three couplings have the values gf {yo )/g'^ = 3.6, g^ (yo )/fi'^ = 
28 and g1 (yo )/fi'^ = 12 ). The blowing up of one or more of the couplings is normally 
interpreted as a sign that the system becomes disordered at that length scale - i.e. at 

Lq ~ a e^o ~ a e^'^/^" ~ a e^'^^ . (3.45) 

Beyond this length scale, the action in Eq. (3.29) is no longer valid and a new action 
(presumably describing massive excitations and, perhaps, some massless excitations) must 
become applicable. Although our calculation cannot be used to obtain any information 
other than the scale Lq at which the perturbative calculations break down, it is interesting 
to speculate on the possible scenarios at longer distances. We can think of three possible 
scenarios. 

(a) All three modes ct, (3 and 7 become massive at Lq , so that there is a gap in the 
spectrum A ~ 1/i^o • 

(b) Only the 7-mode becomes massive at Lq . The a and P modes remain massless until 
a longer distance Li after which they too become massive. 

(c) Only the a and P modes become massive at Lq , but the 7-mode becomes massive only 
at a longer distance Li . 

Clearly, much more analysis is needed to decide between the three possibilities [ 16 ]. 
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Sec. (D): A Large-N Approximation 

A non-perturbative method of studying the mass generation in any field theory is 
to go to the large-A^ limit of an appropriate generalisation of the model. For the 5"^ 
model studied in Sec. 2, the generalisation to an 5"^ model was obvious. Here, however, 
the obvious generalisation of ^ G SO{3) to R E SO{N) runs into problems. To study 
the 5*0 (A^) -valued field theory , it is necessary to impose the constraints using Lagrange 
multiplier fields. But this increases the number of degrees of freedom from A^(A^ — l)/2 
to N"^ , which does not match asymptotically in the large-A^ limit. 

To find a suitable generalisation, let us write the Euclidean SO {3) -valued field theory 
as 

Ce ^ ^[ii + $2 - ( <?i ■ <^2 )' + 0i" + $2 ' ] (3.46) 

where we have eliminated (j)^ = (f)^ + (f)2 — '^{ (pi ■ (p2 )'^ , but we still have to 
impose the constraints (pi^ = (p^ = 1 and (pi ■ (t)2 = 0. Note that Ce still exhibits the 
SO{3) X 5*0(2) symmetry. We can now generalise the three dimensional vector fields (pi 
and (p2 to A^-dimensional vector fields and write the large- A^ Lagrangian as 

Ce ^^[^i^ + $2^ - {$1 -$2? + <?/ ' + 4" ] ■ (3.47) 

This Lagrangian has an SO{N) x 5*0(2) symmetry and since it still satisfies the same three 
constraints, it has 2N — 3 degrees of freedom. We introduce three Lagrange multiplier 
fields to impose the constraints and a fourth field to reduce the quartic term in (3.47) to 
a quadratic form. We then obtain 

^E =^[$1 + $1 + 4>{' + li" + ^Ai (^f - 1) + i\2 {ll - 1) 

V _ _ (3.48) 

-F 2zA3 (p\- (p2 + A4 -h A4 {(p\- (p\ - (p2- (p\)\ ■ 

(By integrating over the \i fields in the path integral, we get back the Lagrangian in 
Eq. (3.47) ). Now, by integrating out the (pi fields, we compute the effective action as a 
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function of A, as 



where 



-S(Xi ) 



'D(f)iD(f)2 exp 

N 



(fx Ce{(I)1,4>2,K 



•^ i=l 1=1 -^ ^ ' 



(3.49) 



M, 



>ii 



i-k) </>2^ (k) ) 



k"^ + iXi iXs + iX^ ko 



1X4 ko 



k^ + zAj 



(l^uik) 
4>2i{-k) 



(3.50) 



so that 



^(A, 



N r (Pk 

+ y J J2^ 



N 



. Ai . A2 

9^ 9^ 



ln{e + '-{Xi + X2) + I \j ^^' ^^' ^' + (A3 + A4fco)2 } 



+ In { A;2 + - (Ai + X, 



iM ^^' /^^' + {Xs+x^kor } 



(3.51) 



Because of the factor of A^ in S{Xi ), in the large- A^ hmit, the integral over A^ in the path 
integral is dominated by the saddle-point of S{Xi ) given by dS/dXi = 0. This is found 
to be at Ai = A2 = —i rn? and A3 = A4 =0 where the saddle-point equation gives 



1 



(fk 



(27r)2 (A;2 + ^2 



Hence 



m = A exp (- 27r/^2 ) = A exp (-2.57 S) 



(3.52) 



(3.53) 



where A is the ultra-violet cut-off a~^ . We can verify that this solution is actually a 
minimum by computing Dij = d^S /dXi dXj at the saddle-point and showing that all its 
eigenvalues are positive. In fact, 



/1/4 



rem 



D^J 






V 



\ 

l/Anm'^ 

l/27rm2 

1/g^/ 



(3.54) 
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at the saddle-point. Substituting the values of A^ at the saddle-point back in Eq. (3.48), 
we find that 

Ce ^ j^[(t)i + <p2 + <Pi^ + <P2^ + rn" {(P,^ + <p2^) - 2m2 ], (3.55) 

so that all 2N fields have become massive. 

Thus, the large- A^ approximation confirms our earlier result that the field theory of 
the large-5' limit of the Majumdar-Ghosh model actually describes a disordered phase with 
a mass gap m ~ e~^^ ^ where K is some number of order one. 
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4. Discussion and Outlook 

Our results in the previous section clearly indicate that the Majumdar-Ghosh model 
is exponentially disordered beyond some distance scale and has a mass gap for all large 
spins S. Thus, for all half-integer spins, by applying the Lieb-Schultz-Mattis theorem 
[ 17 ], we see that the ground state has to be doubly degenerate, with each state breaking 
parity (defined as refiection about a site). For 5" = 1/2, these ground states are just 
the dimerised ground states mentioned at the beginning of Sec. 3. For higher half-integer 
spins, the ground states are not easy to visualise, although the theorem is known to hold. 
However, for integer spins, the degeneracy or non-degeneracy of the ground state is still 
an open question. 

Can the procedure described in Sees. 2 and 3 be used to study general spiral phases? 
The derivation of the field theory presented in Sees. 2 and 3 crucially depended upon 
grouping together all the spins in one period and then identifying 'large' and 'small' vari- 
ables. These became the non-linear fields and their canonical momenta in the long-distance 
field theory. This procedure clearly fails for the general spiral phase for which the classical 
ground state is aperiodic. 
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Figure Captions 

1. Classical ground state of the Majumdar-Ghosh model. 

2. Interaction vertices between Rq and rj. The dotted lines denote the Rq fields and 
the solid lines denote the rj fields. 

3. One-loop effective action for ^q obtained by integrating over the rj fields. 
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